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Abstract 

We study the production and decay of heavy selectrons in the e~7 mode of a hnear 
collider of the next generation. The standard model backgrounds can be substan- 
tially reduced by appropriate kinematical cuts. As a consequence, selectrons far 
heavier than the kinematical threshold for pair production are shown to be easily 
discoverable for large portions of the supersymmetry parameter space. We also de- 
scribe a model-independent kinematical measurement of the mass of the lightest 
neutralino. 



1 Introduction 



One of the most promising ideas that take us beyond the standard model is 
supersymmetry However, while it cures many of the ills that plague the 
standard model, it has its own attendant problems, not the least of which 
is the fact that, till date, we have seen no evidence for its existence. This 
naturally sets a minimum scale for the breaking of this symmetry. A second 
problem deals with the proliferation of parameters in such a theory, especially 
in the context of the soft supersymmetry breaking terms that must be intro- 
duced. Though several attempts have been made to constrain the parameter 
space by demanding consistency with the experimentally measured rates for 
low-energy processes, such bounds are most often easily evaded if the mass 
splittings between supersymmetric particles of a given kind are comparatively 
small. In addition, such constraints have the drawback of being of an indirect 
nature and hence dependent on cancellations (or lack thereof) between the 
contributions from various diagrams with different virtual exchanges. 

A complementary (or even better) method would be to look for an actual 
production of the superpartners in high energy collisions. Of course, the direct 
production of new particles is only possible provided the latter are sufficiently 
light. Typically, a collider facility requires a center of mass energy comfortably 
in excess of twice the mass of the particles to be discovered to permit their 
pair-production. Consequently, the tightest bounds on the weakly interacting 
supersymmetric sector are at present given by the LEP experiments and come 
close to mz/1. 

In this article, we aim to perform an analysis in the context of the linear 
colliders which are currently being planned, such as CLIC, JLC, NLC, TESLA, 
VLEPP, etc.. To be specific, we concentrate on the "canonical" design with a 
center of mass energy = 500 GeV, an electron beam polarization of 90% 
and an integrated luminosity of 10 fb~^. Though most often these machines 
are thought of in terms of e^e~ facilities, they can also be transformed to 
function in the e~e~, 77 and e~7 modes. The last-mentioned is the one which 
is most suitable for our purpose as it lends itself to single selectron production, 
circumventing in this way the kinematic bound that pair production entails 
in e^e~ , e~e~ and 77 collisions. 

The production of a single selectron in association with a neutralino was ana- 
lyzed in Refs We improve these studies here by taking fully into account 
the gaugino mixings, considering polarized electron and photon beams and 
performing a detailed analysis of the differential cross sections of the signal 
and its standard model backgrounds. The latter allows us to devise very effi- 
cient kinematical cuts which dramatically enhance the signal to background 
ratio. As a result, the supersymmetry parameter space can be explored much 



1 



deeper than previously expected [Q]. 



In the next Section we describe the production and decay of a selectron in e~7 
scattering, as well as the dominant standard model backgrounds. In Section 3 
we explain how photon beams can be obtained at a linear collider. We study 
in the following Section the kinematical characteristics of the event distribu- 
tions for the signal and background reactions. Armed with this knowlegde, 
we analyze in Section 5 which regions of the supersymmetry parameter space 
can be explored with the machine described above. Finally, we summarize our 
results in the Conclusion. 



2 Signal and Backgrounds 



We restrict ourselves to the minimal supersymmetric standard model with 
unbroken i?-parity. In most realistic scenarios, the lightest supersymmetric 
particle is the lightest neutralino Xi, which is stable and hence escapes detec- 
tion. We concentrate here on the reaction 

e"7 — > e'Xi , (1) 

and the subsequent selectron decay 

e- e~x\ ■ (2) 



The branching ratio for this decay can be as high as nearly 100%, but it can 
also be substantially lower if other two body decay modes are kinematically 
allowed: 

e--^e-x° (2 = 2,3,4), (3) 
e--^iyeX7 (^=1,2), (4) 

where the Xi and x7 a-^'^ any of the neutralinos or charginos in the theory. 

Obviously both the cross section (|l|) and the partial decay rate (||) depend 
on the details of the supersymmetric extension, particularly the masses and 
the mixing angles in the chargino and neutralino sectors. These receive con- 
tributions not only from terms originating in the naive supersymmetrization 
of the standard model lagrangian, but also from soft supersymmetry breaking 
terms. The chargino mass matrix, in the {W~ h^) basis is given by IQ 
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whereas the neutrahno mass matrix, in the (7 



^ h9,) basis, is 



Mn 



/ Micl + M2sl 
(Ml - M2)c^s, 


V 



(Ml - M2)c^s 
Mi4 + Msg 
mzsp 
-mzcg 



2 
w 





mzsp 




\ 

-mzC/3 

/ 



(6) 



Here Mi^2 are the soft U{1)y and SU{2)l supersymmetry breaking masses 
which parametrize direct gaugino mass terms, whereas /i is the higgsino mixing 
mass. Furthermore, Sw = y^l — = sin^^; and = ^1 — c| = sin^/3 where 
9^, is the weak mixing angle and tan/3 = Vi/v2 is the ratio of the vacuum 
expectation values of the two Higgs fields. We assume the GUT relation 

Ml = ^Ms tan^ 9^ , (7) 



thus leaving us with three parameters namely /j,, M2, tan f3. The exact eigen- 
system has been computed in Ref. We confirm these results and use them 
in our calculations. 

One further question remains unsettled, though. It concerns possible mixings 
between the various charged lepton states. However, both constraints from low 
energy experiments and theoretical studies suggest that such mixings 
are very small. We shall thus assume these to be absent. Also, for simplicity, 
we regard and e]^ to be degenerate, with mass rrie- 

The signal we focus on is thus a final state comprising just a single electron 
associated with missing energy and momentum: 

e"7 — ^ e~Xi — ^ e"x?X? • (8) 



A considerable background exists, though, on account of the standard model 
resonant processes 

e~7 — > e~ — > e^UiUi (9) 



and 



e 7 



w- 



e Up Up 



(10) 



We neglect non-resonant processes which also contribute to the same observ- 
able signature 

e~ 7 — > + Ui + Ui . (11) 
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where z is a generation index. This is an excellent approximation which 
provides a clear picture of the kinematic distribution of the final state. We 
can safely use the narrow width approximation to compute the total and 
differential cross sections of the processes (||^), since the decaying particle is 
either a scalar or it decays invisibly. The subprocess cross sections for these 
reactions have been first calculated by Renard They were also calculated 
for particular helicity combinations in Ref. [^. The expressions in Ref. ^ 
carry some minor errors. 

For the process (JIDD, though, the spin-angle correlations have to be taken 
carefully into account, because the final state electron originates from the 
decay of a vector boson. To do this, we have treated this reaction as a resonant 
2 — i> 3 process. The analytic expressions for the differential cross sections are 
too long to be displayed here. 

Before embarking on further calculations, it is worth noting that the back- 
ground process ([T0|) can of course be very simply reduced with polarized elec- 
tron beams. Unless stated otherwise, we shall consider in the following a 90% 
right-polarized electron beam. 



3 Photon Beams 



High energy photon beams can be obtained by back-scattering a laser ray off a 
high energy electron beam [|TI|. The result of this Compton scattering is that 
the electrons are deflected and dumped while the photons are boosted into a 
hard collimated beam. If the laser is sufficiently intense, all electrons have a 
chance to interact, so that there is no loss in luminosity. 

Consider laser photons of energy Ei {0{1) eV) and circular polarization Pi 
scattering with electrons of energy Ei, ((9(10^^) eV) and longitudinal polariza- 
tion Pf,. (In the next sections we will also need to define the polarization of 
the the other electron beam, which is not Compton-converted.) If the angle be- 
tween the laser and electron beams be 9bi, the scattering may be parametrized 
in terms of 



^ = — ^cos^f' (12) 



where rrie is the mass of the electron. For large z, multiple scattering ( e.g. , 
pair creation in conjunction with a laser photon) becomes important [|12| and 
this results in a depletion of the high energy end of the spectrum. To prevent 
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this, one needs to ensure that 



z < 2(1 + ^2) . 



(13) 



In the following we shall assume the equality in the above expression. The 
resultant photon normalized energy spectrum n[x) and polarization P^{x) are 
depicted in Figs |l] as functions of the fraction x = E^/Ehof the electron energy 



carried by the photon. The underlying analytic expressions are |TT 



dn{x) 
dx 



if 1 4x 4a;2 

AT I ^ ~ ^ ~ z(l -x)^ zHl - xY 



+ PbPi- 



x{2 — x) 
1 — X 



2x 



z{l — x) 



PiC{2 - 2x + x^) + Pbx{l + C 



(1 -x){2-2x + x2) - Ax{z - zx- x)/z'^ - PbPiCx{2 - x) 



(14) 



(15) 



where 



< X < 



z + l 



(16) 



^. 2^ + 18^2^24^ + 8 / 4 8\ ,^ 
^= 2ziz + lY +(^"-^)^^(^ + -^ 



2 - 



{z + lf 
C = l-x(l + l/z) . 



1 



\n{l + z) 



(17) 
(18) 



By design, the energy of the photons can never reach the beam energy since 
their spectrum is limited by Eq. (|13|,|T6|). Furthermore, in practice the low 
energy tail of the photon spectrum cannot participate in any reaction either. 
Indeed, as is displayed in the last of Figs |I], there is a one to one relationship 
between the energy of the back-scattered photons and their angle with respect 
to the direction of the initial electron: harder photons are emitted at smaller 
angles whereas softer photons are emitted at larger angles. For small deflection 
angles from the beam direction we have 
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Clearly, the photons will be distributed according to an effective spectrum, 
which effectively throws out the low energy photons, these being produced at 
too wide an angle to contribute significantly to any reaction. The exact profile 
of this effective spectrum depends somewhat on the shape of the electron 
beam. In the absence of a detailed (and machine specific) study of this effect, 
we approximate this effect by a sharp cut. The position of this cut depends 
of course crucially on the conversion distance, i.e. the distance between the 
interaction point and the point where the laser photons are back-scattered. 
For example, if we assume a conversion distance of 5 cm and an interaction 
spot of 500 nm diameter, then photons scattered at more than 5 /ird are lost. 
The lowest energy photons would then have be read from 

the last of Figs |I[ In the following we use (c/. Eq. (p!3|)) 

Xrain = -5 w = ^ _^ ~ 0.8284 . (20) 

We have checked that our final results are not very sensitive to the choice of 
the first of these two machine parameters. Any cross section is obtained by 
convoluting the fixed energy scattering cross section with the photon spectrum 

my- 



da{s) = J dx ^-da{xs) 6 ix 




(21) 



where the sum runs over the masses of the final state particles. 



4 Kinematical Cuts 



In the absence of linear polarization, the final state electrons in either the 
signal or the background processes are characterized by only two indepen- 
dent kinematical variables, say the cosine of their polar angle (cos6'e) and 
their energy (-Eg)- As it turns out, the distribution of events in this plane is 
very different for the three processes (^p!0|). We now proceed to study the 
boundaries of these distributions in order of increasing complexity. 



4.1 e-j^e-Z^ 



Were the initial state photons monochromatic {E^ = the final state 

electrons would be distributed in the (cos6'e, i^^e) plane along the line of con- 
stant energy Ee = {s — m%)/2^. Since the energy of the initial state photon 
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is actually spread between XminV^/2 < < Xma.x\/s/2 (]20|), the final state 
electrons are distributed in phase space region 

n _ (1 + a^)v^^e -xs + ml 

(1 - x)^/sEe 



The two boundary curves are displayed in Figs |^ and |^. They confine all the 
events, and will be used as kinematical cuts in the subsequent numerical 
analysis to entirely eliminate the background (P). 

4.2 e-7 e'Xi e'XiX^ 



The electrons emerging from the decay of the selectron are distributed within 
the phase space region 

2 2 
m'l — m| 

' = 2[E,-k,cosie,-9,)] ' ^^^^ 



where kg = y E? — m| is the momentum of the selectrons, whose energy 
spreads within the region 

_ (1 + x)y/sEg -xs-ml + m\ 
(1 - x)^/skg 



Since the (heavy) selectrons are produced rather isotropically 0, the cosine 
on the left hand side of Eq. can take any value between —1 and 1. As a 
consequence, the highest energy of the selectrons is given by 



^5 



+ (1 - Xmax) l/ ( a;maxS + TTlj - m|) - 4XmaxSm? 



(25) 



and the energy of the decay electrons is confined within 



2 2 2 2 

m| — mi mi — mi , , 

— ^ <Ee< — — ^ — - . 26 

e 



2 (Ef ^ + fcf''^) - 2 (E?'^^ - fc, 



The two boundary lines are displayed in Figs and |^ for a selectron of mass 
mg = 250 GeV and the particular set of supersymmetry parameters /i = 500 
GeV, M2 = 400 GeV and tan = 4. For this choice, the mass of the lightest 
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neutralino is m^o = 193 GeV. The available phase space is thus small because 
the sum of the produced masses rrie + m^o = 443 GeV is so close to the 
available center of mass energy -^/Xmax-s = 455 GeV. 

Note that the values of the two energy boundaries ( ^61) depend solely on the 
center of mass energy, the mass of the selectron and the mass of the neutralino. 
Therefore, a measurement of these two limiting energies provides a direct 
kinematical determination of the mass of the neutralino. A similar model 
independent measurement of the neutralino mass can otherwise be performed 
only in a polarized e~e~ experiment, if the selectron is light enough to be 
pair-produced |T^. 



4-3 e 'J W e Ue^e 



Here the situation is more subtle. At first glance one might expect the decay 
electrons of the W~ to be subjected to similar kinematical bounds as the 
decay electrons of the selectrons. Indeed, proceeding with the replacements 



mg ^ my/ and ^ niy = in Eqs (p3| , p6D , the energy of the W decay 
electrons is in principle only bound by 

0<^,<f. (27) 



This range covers almost the entire phase space. However, some improvement 
can be achieved by identifying the dense regions in phase space. Since the W~ 
production cross section is very much peaked backwards with respect to the 
initial electron beam [Q, it is a good approximation to set co^Ow = —1 in 
Eqs ( p3| , |2^ ) . Hence, the most probable energy of the W~ is 

2 2 

_ X s + ruw < T < T ('?R\ 



and the bulk of the electrons is produced within the phase space region 

< . (29) 

2[Ew + kw cosOe) 



Since the W~ production cross section increases monotonically with the center 
of mass energy more events occur for larger values of x. In Figs ^ and 

^, we have thus displayed the curve corresponding to x = x^ax in Eqs (^ , ^91) . 
Clearly, a large portion of the W~ events is confined below this curve. It will 
therefore be used as a kinematical cut alongwith Eq. (p2D in the subsequent 
numerical analysis. 
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It should be borne in mind that events which do not satisfy this bound are also 
kinematically allowed. However, the corresponding rate is strongly suppressd 
by the dynamics of the W~ production mechanism. One effect of the dynamics 
can be seen by comparing the density of events plotted on Figs |^ and ^. The 
only difference consists of different initial state polarizations. Although the 
choice of Fig ^ confines the W~ events better within the bounds ( ^^9]) , far 
fewer selectrons are produced and the signal to background ratio is worse. 
This is why we concentrate in the next section on the choice of polarizations 
used in Fig El. 



5 Results 



In addition to the kinematical cuts (p2| , p9D discussed above and depicted in 



Figs I and ^ the signal to background ratio can be further enhanced by divid- 
ing the (cos^e; -E'e) phase space into bins (of equal size, here). The number 
of events in each bin is then compared with the standard model expectation. 
This procedure takes automatically into account the information contained in 
Eq. (^61). The significance of the deviation is given by the t^st 

x' = E{^;^) . (30) 



where n is the number of events expected for the standard model, n'^^^ = 
n{fi, M2, tan (3, mg) is the corresponding observed number of events (if super- 
symmetry is to explain the deviation) and the error An is the quadratic com- 
bination of statistical and systematic errors on the observed number of events 

An = ^n + (en)2 . (31) 



The relative systematic error e is essentially due to the luminosity measure- 
ment (the uncertainty on the electrons' energies and angles is negligible) and 
is set to 1% in the following. 

Clearly, if less than five events are contained in a bin, the probabilistic inter- 
pretation of the fsst becomes unreliable. Indeed, with so few events the 
underlying Poisson distribution does not resemble enough a gaussian shape to 
warrant the sum (^) to be distributed according to a x^- For this reason, we 
ignore alltogether any bin with less than five events. 

Instead of attempting to explore at one go the four- dimensional parameter 
space viz. (/x, M2, tan/3, mg), we choose to present our results in the form of 
2-dimensional contour plots in the (/x, M2) plane for different values of mg 
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and tan (3. We also exhibit the dependence on the polarization of the electron 
beams Pb and Pg, and on the number of bins A^. To obtain the explorable 
regions at the 95% confidence level we set = 6 in Eq. pO|). 

Unless stated otherwise we use in the following 

me = 250 GeV tan /? = 4 , (32) 

Pi = -1 = 100% left Pb = .9 = 90% right Pe = -9 = 90% right , (33) 

iV = 3 X 3 . (34) 



Since we consider a machine operating with a e~^e~ center of mass energy 
V^max ~ GeV, selectrons with a mass up to 250 GeV can in principle be 
pair-produced and detected in the e~^e~ or e~e~ colliding modes. We there- 
fore display in Fig. ^ observability contours for rrig > 250 GeV, i. e. beyond 
the kinematical limit of the parent coUider. As expected, the explorable area 
shrinks with increasing selectron mass. The sharp drop for small /i can easily 
be understood as for such values, the lightest supersymmetric particle is pri- 
marily a higgsino and (almost) does not couple to the electron. This part of 
the parameter space can be better explored using the e"'"e~ mode to observe 
chargino pair-production. The regions of parameter space located between the 
dotted lines are those which could be explored this way at the same collider 
(or have already been clearly excluded by the LEP experiments) assuming 
charginos can be observed up to the very threshold for their pair-production. 
If this limit can indeed be reached in practice, the e~7 mode has a chance of 
(izscoi'ermg' supersymmetry where the e~^e~ mode cannot, only if the selectron 
mass is comprised within 250 < mg < 325 GeV, and this in a marginal portion 
of the supersymmetry parameter space. Still, even if supersymmetry has been 
discovered before, for a 500 GeV machine, the e~7 collider operating mode 
is the only one susceptible of discovering a selectron heavier than 250 GeV. 
Even a selectron as heavy as 400 GeV can be observed and studied. 

In Fig. ^ we display the tan j3 dependence. As is easily evinced, the contours 
tend to become more symmetric for progressively larger values of tan (3. The 
reason can be understood by examining the characteristic equation of the 
matrix (^. In the limit where tan/? = oo there remains only a /i^ dependence. 
Similarly, for 3> mz the dependence on tan/3 should vanish. This fact is 
indeed reflected by the convergence of the contours for large values of 

Some portion of the W background (PUD being irreducible, it is important 
to reduce it as much as possible by polarizing the initial electron beam. As 
can be inferred from Fig. |^ the incidence of a poor polarization is particularly 
dramatic in the region 50 GeV ^/z^ 200 GeV, where the supersymmetric signal 
(j^) is small on account of the lightest supersymmetric particle being primarily 
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a higgsino. Nevertheless, improving the polarization beyond 90% only yields 
marginal improvements. This is because even though the background is further 
decreased, the signal remains at the hmit of observability. The erratic features 
of the curve at 95% polarization reflect our rejection of bins containing less 
than five events. 

Finally, in Fig. ^ we turn to the dependence of our resuhs on the binning, 
more specifically on the number of bins. Since coarser binning loses infor- 
mation about the differential distribution, the improvement in the bounds 
with better energy and angular resolution is not unexpected^. It should be 
noted, though, that increasing the number of bins results only in a modest 
increase in sensitivity. This is a direct consequence of the efficient kinematic 
cuts Eqs OH). 



The variation of this improvement with /i can again be traced to the neutralino 
mass matrix (^). For small the lightest neutralino has a large higgsino com- 
ponent. As a result, its coupling to the electron and hence both the production 
cross section (eqn. |^) and the partial decay width (eqn. are suppressed. This 
leads to a comparatively smaller signal to noise ratio. To offset this loss, ad- 
ditional information as obtained from binning is useful. For large values of 
though, the higgsinos tend to decouple and the lightest neutralino is primarily 
a gaugino. The signal to background ratio on imposition of the kinematic cuts 
is sufficiently large to render binning almost inconsequential. 



6 Conclusions 



We have studied the production and decay of heavy selectrons at a linear 
collider of the next generation operated in its e~7 mode in the context of the 
minimal supersymmetric standard model. While the other modes (e^e~, e~e~ 
and 77) are limited by the kinematical limit of rrie < -\/s/2 = 250 GeV, the 
e~7 option can discover selectrons which are much heavier, up to 400 GeV. 
The standard model background can be controlled by a judicious choice of 
beam polarizations and kinematical cuts. 

The phase space distribution of the signal final state electrons stands out 
significantly from the background. This allows to infer in a model-independent 
way the mass of the selectron as well as that of the invisible lightest neutralino, 
the lightest supersymmetric particle. 



This improvement is not a monotonia function of bin cardinality, as very fine 
binning would leave too few events in each, thus disqualifying them from contribut- 
ing to the function. As a matter of fact, a 4 x 4 binning already degrades the 
resolution! 
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X 

Fig. 1. Back-scattered photon energy spectrum, polarization and polar angle as 
functions oi x = E^/Eh for three different combinations of beam polarizations. 
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COS 6e 

Fig. 2. Scatter plot showing the angle-energy distributions for the final state elec- 
trons emanating from the processes (p|,|g| jlO|) . Assuming an integrated luminos- 
ity of 10 fb~^ each points corresponds to one event. The curves show the ex- 
act {Z^,e~) or approximate {W~) kinematically allowed ranges. The laser pho- 
tons are 100% left polarized whereas both electron beams are 90% right polarized 

{Pl = -1 Pb = +.9 Pe = +.9). 
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Fig. 3. Same as Fig. |2|, except for the laser and Compton-converted electron beams 
which have now opposite polarizations {Pi = +1 = — .9 Pe = +.9). 
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Fig. 4. Contours of = 6 for selectron masses mg = 250, 300, 350, 400 GeV (from 
upper to lower curves). All other parameters are specified in Eqs (^^34). The areas 
below the plain curves can be explored with 95% confidence. The dotted curves 
delimit the region already excluded by LEP I (lower curves) and the region below 
which charginos can be pair-produced at the same linear collider run in the e~^e~ 
mode (upper curves). 
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Fig. 5. Contours of = 6 for tan /3 = 1,4, 50. All other parameters are specified in 



Eqs (32-34). The areas below the curves can be explored with 95% confidence. 
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/i [GeV] 

Fig. 6. Contours of = 6 for electron beam polarizations Pb = Pe = 0, 50, 90, 95% 
(from lower to upper curves). All other parameters are specified in Eqs (|32|- 34 ) . 
The areas below the curves can be explored with 95% confidence. The dotted curves 
depict the kinematical limit of the supersymmetric reaction (|8[). 
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Fig. 7. Contours of = 6 for a number of bins N = 1x1, 2x2, 3x3 (from 
lower to upper curves). All other parameters are specified in Eqs (^^34). The areas 
below the curves can be explored with 95% confidence. The dotted curves depict 
the kinematical limit of the supersymmetric reaction dq) . 
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